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1 Introduction 


This document is meant as a pedagogical introduction to the modern language used to talk 
about quantum theories, especially in the field of quantum information. It assumes that the 
reader has taken a first traditional course on quantum mechanics, and is familiar with the 
concept of Hilbert space and elementary linear algebra. 

As in the popular textbook on quantum information by Nielsen and Chuang [I], we introduce 
the generalised concept of states (density matrices), observables (POVMs) and transforma¬ 
tions (channels), but we also go further and characterise these structures from an algebraic 
standpoint, which provides many useful technical tools, and clarity as to their generality. This 
approach also makes it manifest that quantum theory is a direct generalisation of probability 
theory, and provides a unifying formalism for both fields. 

Although this algebraic approach dates back, in parts, to John von Neumann, we are not aware 
of any presentation which focuses on finite-dimensional systems. This simplification allows us to 
have a self-contained presentation which avoids many of the technicalities inherent to the most 
general C*-algebraic approach, while being perfectly appropriate for the quantum information 
literature. 


2 States and effects 

2.1 Basic Quantum Mechanics 


Let us start by reviewing the standard 1930’s formulation of quantum mechanics, which is still 
used in many textbooks. The theory is defined around a Hilbert space 7/, whose normalised 
vectors | ip) G V. represent states, or “wave functions”. The observables of the theory are 
represented by the set of self-adjoint operators = A acting on %. 

The operational meaning of the observable has to be “unpacked” using the spectral theorem. 
For instance, if "H has a finite dimension d , it states that there exists an orthonormal basis 
{|i)}f =1 , i.e., (i\j) = Sij, such that A\i) = ai\i), where ^ 6 R are the eigenvalues of A. 

In the simple case where the eigenvalues a* are all distinct, the operational content of the 
formalism is summarised in the following statement: when measuring the observable A on a 
quantum system in state |"0), the probability of obtaining the ith outcome (where the observable 
takes the value a*) is pi = (i\if;) = |(i|^)| J . 

More generally, if the operator A has degenerate eigenvalues, it is convenient to write A\i,k) = 
ai\i, k ), where k labels the degeneracy of eigenvalue a^, and the a, are all distinct. The proba¬ 
bility to obtain the ith outcome is then determined by the famous 
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2 States and effects 


Definition 2.1 (Born Rule). 

If the system is in state | ip), a measurement of the observable represented by the self-adjoint 
operator A returns the value a.* with probability Pi = YPk |(b ^\ip)\ 2 . 


Normally, one goes on to describe the projection postulate, but we will see that it can actually 
be derived from the born rule, once properly generalised. 


2.2 Positive operators 


In the first part of the lecture we are going to generalize this description of quantum mechanics, 
and in some sense simplify it. For this purpose, it will be useful to understand the concept of 
positivity of operators. 


Definition 2.2. 

An operator A is called positive, denoted as A > 0, if for all \ip) e %, (ip\A\ip) > 0. We also 
write A>B if A-B>0. 


In finite dimension, this definition is equivalent to A being self-adjoint and having non-negative 
eigenvalues, that is ,/R — A, Vi : a* > 0. It is straightforward to see that such an operator 
implies the defining properties. Let {|z)} denote the eigenbasis of A: 

{ip\A\ip) = Ycl c j{A a j\j) = N 2( b > 0. 

i,j i 

For the converse, consider A with (ip\A\ip) > 0 for all \ip) G B. To see that this A is self-adjoint, 
observe that in general, 

{fp + <f>\A\ip + 4>) = {ip\A\ip) + {<j>\A\4>) + (ip\A\(fi) + {<p\A\ip). (2.1) 

Since this is real by assumption, and also the first two terms on the right side of the equality 
are real, we must have 

hP\A\<t>) + {<f)\A\^}) = (pp\A\4>) + (0|A|^) 

for all states ip and 0, which implies Im (ip\A\(p) = —Im(0|A|^). Replacing (p by i<p, we also 
obtain Re (ip\A\<p) = Re (0|A| ip). It follows that (^|A|0) = (cp\A\ip), and hence A is self-adjoint. 
The positivity of the eigenvalues follows from the fact that eigenvalues can be considered as 
expectation values of eigenvectors, i.e. a* = (i\A\i) > 0. 

Another important point to notice is that for any operator B, we have B^B > 0, because for 
all ip, (ip\B^B\ip) = ||R|r/?) || 2 > 0. It is also true that any positive operator A can be written in 
this form by using, for instance, B = y/A, defined as having the same eigenvectors as A, but 
the square root of the eigenvectors of A. 
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Also, we will use the fact that the binary relation A > B if B — A > 0 defines a partial order 
on the set of all operators. This order is partial because not every pair of operators A, B can 
be compared. For instance, consider the matrices 


A 



and B 



It is easy to see that neither A — B nor B — A are positive, since they both have a negative 
eigenvalue. 


An important set of positive operators are the projectors: 


Definition 2.3. 

An operator P is a projector, if it satisfies P'P = P or, equivalently, P = P' and P 2 = P. 


It follows that the eigenvalues of a projector P are either 0 or 1. Indeed, if P\ %jn) = p\ijj) for 
some eigenvalue p, then also p\ip) = P|-0) = P 2 \ ijj) = p 2 \ip), it follows that p 2 = p, and hence 
p — 0 or p — 1. In particular, this implies that P > 0. 

Such an operator P projects all vectors orthogonally on the subspace PPL = (P|^) : |^) £ 
PL} = {|"0) G PL : P\ip) = |^)}. Given any normalised state \(j>), we can define the rank-one 
operator 

p* = WA 

as a projector on the subspace spanned by \(f>). It maps any vector |^) to P^\i]j) = (0|^)|0). 

This notation allows us to write any self-adjoint operator A (in finite dimensions) as a sum of 
projectors onto its eigenspaces. Let | i, k) be the eigenvectors of A with distinct eigenvalues a*, 
i.e., such that A\i,k) = a.i\i,k ). We define the spectral projectors 

Pi := J2 \i,k)(i,k\ 

k 

which projects on the eigenspace P t PL associated with the eigenvalue a*. (It is easy to check 
that a sum of projectors is also a projector). This allows us to write A in terms of its spectral 
decomposition: 

A = pPf 

i 

The completeness of the eigenvectors imply that P; = 1, although one of the eigenvalue a, 
may be equal to zero, eliminating one of the terms from the sum above. 

Observe that the probability of obtaining the value a* in a measurement of A on state \ij}) can 
now be written simply as 

Pi = Y, l(*>*#)| 2 = (^1 PM>)- 

k 

In particular, it depends solely on the spectral projector P, t and on the state | ijj). 

Note that this simple expression holds only provided that we make sure that the eigenvalues a* 
are distinct and hence that the projectors P* are of maximum rank. Moreover, this makes the 
spectral decomposition of a self-adjoint operator unique. 
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2 States and effects 


Often, one is not interested in the individual outcomes probabilities when measuring an ob¬ 
servable, but just in the average of the measured value: 


Definition 2.4. 

The expectation value (A) of an observable A with spectral decomposition A = Xu a iP% 
measured on a quantum system in state \ij}) is given by 

(A) = Y ciiPi = Y, a i(' l P\P i \' t P) = \A\^) (2.2) 


It is worth noting that the probabilities p % themselves are expectations values of the spectral 
projectors: 

Pi = (Pi)- 

Therefore, any prediction of quantum theory is given by the expectation value of some self- 
adjoint operator. 


2.3 Generalized States 


In what follows, we introduce two settings which show that the concept of a state needs to 
be generalized from vectors in a Hilbert space to certain operators with special properties. 
We start by giving two situations motivating this generalisation, and then give an abstract 
definition which encompasses both examples. 


2.3.1 Ensembles 


Consider the case where the observer of quantum system is uncertain about what its exact 
state is. A natural way to model this situation is to assign probabilities p t to quantum states 
i/)i according to his belief about the system. This defines the ensemble {(pi, V’*)}?= 1 - 

The expectation value of the observable A must then be the average, in terms of the classical 
probability distribution i H» p^, of the various quantum expectations values 


{A) = YPi^Ml^i)- ( 2 - 3 ) 

i 

This can be rewritten in a more compact form using the trace operator Tr. Recall that the 
trace of a matrix is cyclic, i.e., Tt(AB) = Ty(BA). Of course, the product AB must be a square 
matrix, otherwise the trace is not defined. However, A and B themselves can be rectangular 
matrices. The best way of thinking about a ket \^i) is as a matrix with just a single column, 
whereas a bra |?/?;) is a matrix with just a single row. Considering that the trace of a number 
(i.e. a one-by-one matrix) is that number itself, we obtain 

{ipi\A\ii>i) = Tr(^| = Tr(|*/'i)(t/'ilA 


(2.4) 
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We invite the reader to verify that this is correct. This allows us to rewrite the expectation 
value of A as 

(-4) = \ A \^i) = 52Pi Tr(l^iX^il^) = T r(pA) (2.5) 

i i ^ ' 

where we have defined the operator 


p ■= Hi’ilV’iXV'i 


( 2 . 6 ) 


which is usually called density matrix , or density operator. Given that, as noted at the end of 
Section [272], all quantum predictions take the form of the expectation value of some self-adjoint 


operator, this means that in this scenario, the matrix p is all that we need to know about the 
ensemble {(p,,'0*)}?=! i n order to compute predictions. 


For instance, considering an observable in its spectral decomposition A = Jfj a jPji the proba¬ 
bility Qj of observing the outcome a 3 is also simply 

(Pi) = Tt(pPj)- ( 2 . 7 ) 


It is important to not that many different ensembles give rise to the same density matrix p : 


Theorem 2.1. Two ensembles {( Pi , |Vk))}"=].• {(Qu |0i))}£li are represented by the same density 
operator, i.e, 

J2Pi\A)(H = ( 2 - 8 ) 

i i 

if and only if there exists an unitary matrix Uij such that 

VPi\A) = ( 2 - 9 ) 

3 


Proof. The sufficiency of the condition can be verified by simple substitution and using the 
unitarity of the matrix u t j, namely the fact that JfiUikUij = dkj- For the converse, we follow 
Ref. pp. Let p = PilV^KV^I- First, observe that we can always build an ensemble of mutually 
orthogonal states representing p using its eigenvectors | k) and eigenvalues pp p = JfiPi |*)(*|. 
If we can prove the theorem in the case where one of the ensembles is composed of orthogonal 
states, then we obtained the desired unitary matrix by multiplying the unitary relating the 
vectors of the first ensemble to the orthogonal ones with that relating the orthogonal vectors 
to the vectors of the second ensemble. Let us define \af) := ^/pil'f’i) and | bf) := y/qi\^>i) for 
conciseness. We assume that the states |a)* are all mutually orthogonal. In this case, observe 
that if \if) is orthogonal to all vectors |a)j, then 

o = (V’HV’) = EKWI 2 

i 

which implies that {bfif) = 0 for all i, and hence \i(j) is also orthogonal to all vectors | bf). This 
means that the families | af) and | bf) have the same span. Hence there are complex numbers c t j 
such that 

| b i ) ) ( Cij | (I j ). 

3 
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2 States and effects 


Moreover, we have 


P = E W){a k | = E l 6 *X 6 fcl = E E Cki^kj | ®i) (®fc | • 



From the linear independence of the matrices |a;)(a0, we conclude that J2k c kiC k j = $ij, and 


hence the matrix c tJ is unitary. 


□ 


The uncertainty involved in a given ensemble can be measured by the Shannon entropy S(p) = 
— Y^i 'Pi hip,> For a given density matrix p, this entropy is not defined because it depends on 
the ensemble from which p is constructed. However, it makes sense to consider the ensemble 
that corresponds to a minimal uncertainty. This defines the von Neumann entropy associated 
with a density matrix: 


S(p) = — Tr(p In p) = min S(p), 

P=J2iPi \H)(A\ 


( 2 . 10 ) 


where the minimum is over all possible ensembles {(p ?; ,0 ?; )}” =1 such that p = Y^iPi\' l l , i)(' l l , i\- hi 
fact, one can show that this minimum is reached whenever the states 0j are all orthogonal. In 
this case, the probabilities pt are simply the eigenvalues of p. Therefore, S(p) is the Shannon 
entropy computed from the eigenvalues of p. We will come back to this special diagonalisiug 
ensemble in Section [2.3.31 


2.3.2 Tensor Product and Reduced States 

For this example, we need to introduce another tools in quantum theory which allows us to 
compose several systems into a bigger one, or, alternatively, to consider subsystems of a larger 
system. 

Consider two quantum systems A and B with respective Hilbert spaces 1-La and 1-Lb■ Let 
{|i)u}" = i be an orthonormal basis of Ha and {|j)bKLi an orthonormal basis of Hb respectively. 
We define the larger Hilbert space 


( 2 . 11 ) 


n = a a ®Bb 


of dimension n ■ m, with orthonormal basis — — m. It is convenient to 

write formally 


i,j) = NU® Ij)b. 


( 2 . 12 ) 


Hence, any vector -0 G H can be expanded as 


n m 




(2.13) 


i= 1 j= 1 


We can think of Ha and Hb as two “independent” parts of H. In particular, given any vectors 
£ Ha and 0# G Hb, we can construct the joint vector 


n m 


|0u) ® 10b) = X!Z!(*IVu)O'l0iv> I*)a® \ j ) b - 


i =1 j= 1 


(2.14) 
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Moreover, there is a natural way to map any operator A acting on "Ha, or B acting on H B , to 
an operator on H, respectively written as A ® 1 and 1 ® B. They are defined by their action 
on the basis as follows: 


(A ® l)(|i)A \j)b) = (A| i) A ) ® | j) B (2.15) 

(1 ® B)(\i) A ® | j) B ) = \i) B ® (B\j) B ). (2.16) 

Observe that the labels A and B on the kets are not really necessary, as we always take care 
to preserve the ordering of the tensor factors: system A on the left of the tensor product and 
system B on the right. Hence in the following we would write the above two equations as 

(A ® l)(|i) ® | j)) = (A|i)) ® | j) (2.17) 

(1 B)(\i) ® | j)) = | i) ® ( B\j )). (2.18) 


An essential property of this representation of the operators A and B is that they commute: 


(A ® 1)(1 ® B) = (1 ® B){A ® 1) = A ® B. (2.19) 

We also have the following algebraic properties, 

(A ® B)(\iP) ® |0)) = (A|0)) ® (H|0)) (2.20) 

(A <g> B){A! ® B') = (AA') ® (BB') (2.21) 

a(A® B) = (oiA) ® B = A® (aB) (2.22) 

(A ® H) f = A f ® B f , (2.23) 


for any operators A, A', B , acting on the respective Hilbert spaces, and any a E C. 

The tensor product can be straightforwardly generalised to more than two tensor factors, and 
two non-square matrices, i.e., operators between two different Hilbert spaces. In particular, the 
tensor product of two kets is the same as that of two matrices with only one column (columns 
vectors). We invite the reader to experiment with this concept. In particular, a good exercise 
is to understand the tensor product in terms of its action on matrix components. Here we just 
observe that the tensor product is associative: 

(A® B) ®C = A® (B ®C) = A® B ®C. (2.24) 


Also, a construction that we will often encounter is the tensor product of an operator and ket, 
or a bra, whose action on states is defined as follows: 

(A ® |0))|i) = A\i) ® |0) (2.25) 

(A ® (0|)(|i) ® | j)) = {ip\j)A\i). (2.26) 


We are now in measure to define the concept of reduced state. Suppose that, for some reason, 
we decide to only measure observables on Ha, he., observables represented by operators of 
the form A ® 1, where A is any operator on Ha- Observe that, in particular, if the spectral 
decomposition of A is A = J2i a iPi, then the spectral decomposition of A® 1 is A = Yli a iPi ® 1- 
Indeed, it is easy to check that Pi® t are also projectors. 




2 States and effects 


This implies that all the possible predictions are of the form of the expectation value of a 
self-adjoint operator having the shape X <g) 1, which, on an arbitrary state T G Ha ®Hr- is 

(X <g> 1} = (ip\X ® 1\4>) = Tr ((X 0 1MM) = £Tr ((* ® \j)(j\) iV’XV’l), (2.27) 

where, in the last step, we expanded the identity on system B in terms of a basis with elements 
| j). Observing that 

X ® \j)(j\ = (1 ® |j»X(l ® (j|), 
and using the cyclicity of the trace, we obtain 

(X ® 1) = 5>((i ® \J)) X ( 1 ® (j\MW) 
j ' 

— X] Tr^X(l ® (jj)|'0)('0|(l ® \j))\ 

3 ' ' 

= Tr(A» 

where we defined the operator 

P = E( 1 ® 01)I^)(V’I(1 <g> | j)). 

j 

This operator p, which is also a density matrix , is an operator acting only on Hilbert space Ha- 
Nonetheless, it contains all the information that we will ever need about the state ip G Ba^Hb, 
provided we are restricted to only measure observables of system A. 

It will be useful to define more generally the operation which maps \'p)('p\ to p. 
map that we call the partial trace over system B , written Tr^, which acts on 
follows 

Tr B (Z) = £(l®0j)Z(l<g)|j>). 
j 

It takes an operator Z acting on 1-La ® Hb to an operator acting on Ha alone. 

This is called partial trace because the full trace Tr over Ha <S> Hb is given by 

Tr = Ti '4 Tt b = Tr B Tr^, (2.32) 

where the implicit product used here is simply the composition of maps. 

The above reasoning can be carried likewise if we started with a density matrix pab on system 
Ha ®Hb rather than simply the vector 


It is a linear 
operators as 

(2.31) 


(2.28) 


(2.29) 


(2.30) 


Definition 2.5. 

Let pab be a density matrix on a bipartite system, then 

Pa = Tr b(pab) (2.33) 

is the reduced state on system A. 






2.3 Generalized States 


9 


2.3.3 Density Operator/Matrix 


In the previous two sections, we have seen that it is useful to characterize the state of a system 
through an operator p, such that the expectation value of a self-adjoint operator A can be 
obtained via Tr(pA). Now we take this more general concept of states to its full generality: 


Definition 2.6. 


A density operator (state) is an operator p satisfying 


0 P > 0, 


ii) Tr p = 1. 


The expectation value of the observable represented by the self-adjoint operator A with respect 

to the state p is given by 


(A) = Tr(pA). 

(2.34) 


These two conditions completely capture the concept of a density matrix that emerged in either 

2.3.2| That is, not only do the density matrices 


of the two contexts studied in Section 2.3.1 


or 


emerging in these context always satisfy conditions (i) and (ii), but also any matrices satisfying 
these condition can emerge in both contexts. 


First, observe that these conditions are necessary and sufficient for the following interpretation: 
for any observable represented by a self-adjoint operator A with spectral decomposition A = 
J2i a iPii we want the numbers p* = Tr(pPj) to form a probability distribution, i.e., that p* > 0 
and YliPi — 1- Indeed, this requirement implies in particular that for any vector |-0) we must 
have Tr(p|^)(^>|) — (fi’lpl'fi) > 0, which simply means p > 0, namely Condition (i). Moreover, 
from the fact that YliPi — 1 , YliPi — 1 directly implies Tr(pl) = 1, which is Condition (ii). 
The converse is similarly straightforward. 


Moreover, any such density matrix can be obtained either as an ensemble or as a reduced 
state, which shows that this definition is not taking us away from the accepted framework of 
quantum mechanics. Indeed, suppose p is any matrix satisfying conditions (i) and (ii). Let |z) 
be a complete set of eigenvectors for p with eigenvalues p;, i = 1,... ,n. The two conditions 
imply that Pi > 0 and YliPi — 1, therefore p represents the ensemble {p,, |f)}” =1 via its spectral 
decomposition p = J2iPi\i){i\- This arbitrary density matrix can also represent a reduced state. 
Indeed, consider two copies of the Hilbert space on which it is defined, and the bipartite vector 
If/’) — Y2i -\fPi\i) <H> |*)- One can check that p is obtained by tracing out \'ip)('ip\ on any of the two 
systems. The vector \fi) is generally called a purification of p. 
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2 States and effects 


Definition 2.7. 

We call a state p pure if it cannot be written as an ensemble of two or more distinct states, 
(i.e., it is extremal in the convex set of density matrices). The following propositions are all 
equivalent (on finite-dimensional Hilbert spaces): 

i) p is pure 

ii) S(p) = 0 

iii) p 2 = p 

iv) p = \ip)(ip\ f° r some normalised vector ip. 

Proof. We show (i) (ii) (iii) => (iv) =>■ (i). It directly follows from the extremality 
condition that S(p) = 0. If p = J2iPi\i)(i\ for orthogonal eigenstates |i), then we have S(p) = 
— EiPi hi pi = 0. But each term —p t In p, is positive since 0 < Pi < 1. Hence, for all i we have 
—Pilnpi = 0, which is true either if p* = 0 or Inp* = 0, i.e., Pi = 1. Hence p is a projector: 
p 2 = p. Moreover, since J2iPi — 1 then only a single p* can be nonzero. It follows that p = |i)(i|. 
Moreover, such a state is extremal because \ip)(ip\ = ppi + (1 — p)p 2 implies, by multiplying 
both sides by (ip | on the left and | ip) on the right, that 1 = p(ip\pi\ip) + (1 — p)(ip\p 2 \ip), or 
p(l — (ip\pi\ip)) + (1 — p)(l — (ip\p 2 \ip)) = 0. Since both terms are positive we must have, 
(ip\pi\ip) = 1, which implies that Pi = P 2 = \' l P)(' l P\- □ 

We close this section with some examples: 


2.3.4 Example: Entanglement 


Consider a bipartite system, with Hilbert space PLa ®PLb- Let {|?)a}”=i an( f {|.7) b}]=i de¬ 
note orthonormal bases of PLa and PLb respectively. The pure state |O) = \i)a ® \i)B 

corresponds to the density matrix 

pab = |*)b)(0'U<8 ) 01b) = 1 H N)01a® I*)01b- 

n i,j n i,3 

The reduced state on system A is obtained by the partial trace: 


Pa 


Tr b(pab) = ® ( k \) pab(i ® \k)) 

k 

k,i,j 

m = Wi*x*i = 


1 

n 


(2.35) 

(2.36) 

(2.37) 


The reduced state pa = 1 jn is referred to as the maximally mixed state. Because it is invariant 
under any unitary transformation U, it is also the reduced state of the states of the form 
(U 1)|H), which just amount to a different choice of basis on system A. Observe that, 
although the full state pab = is pure, and hence S'(|O) (Ci|) = 0, i.e., we possess maximal 
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information about it, its “part” pa has maximal entropy: S(pa) = In n, which means that we 
know absolutely nothing about the state of system A. This is completely contrary to classical 
systems, where knowing the state of the whole system implies also complete knowledge of its 
parts. Quantum states having this non-classical property are called entangled. Hence, we 
call any pure state pab = \^){M of the compound system entangled whenever its parts have 
non-zero entropy. Since the parts of Q) have maximal entropy, |fi) is maximally entangled. 


2.3.5 Example: Bloch Sphere 


We characterise the complete set of states one the two dimensional Hilbert space % — C 2 . For 
this purpose it is useful to introduce the Pauli matrices {cr 0 , cy, cr 2 , <v 3 }, which form a basis of 
the complex vector space of two-by-two complex matrices: 


0o 




0 -i \ ( 1 0 \ 

i 0 ) ’ a3 ~~ { 0 -1 ) ' 


(2.38) 


The set of Hcrmitian matrices with trace one is completely parametrised with three real numbers 
r i, r 2 , r 3 as: 

P = Wo + ^(ncri + r 2 cr 2 + r 3 cr 3 ) = ]- ( 1 r . 3 r j U2 V, n,r 2 ,r 3 G M. (2.39) 

2 2 2 \ ri + ir 2 1 - r 3 J 

Moreover, for the operator p to be positive, its eigenvalues Ai, A 2 need to be non-negative. Since 
we already guaranteed that they sum to one, they cannot be both negative. Therefore, the only 
extra condition required is that their product be non-negative, that is, 

detp = AiA 2 = J(1 -r| - {r\ - ( ir 2 ) 2 )) = ^(1 - rl - r 2 2 - r 2 3 ) = ^(1 - ||f1| 2 ) > 0. (2.40) 

From this we can see that the set of possible states (density matrices) corresponds via the 
above parametrisation to the 3-dimensional ball characterised by ||r|| 2 < 1. This is manifestly 
a convex set whose extreme points; the pure states, lie on its boundary: the unit sphere. 


2.4 Generalised propositions 


In the previous section we generalized the notion of state. We will also need to generalize the 
concept of observable. For this purpose, we first consider the most elementary notion of yes/no 
observables, i.e., observables with only two distinct outcomes. Following the tradition of formal 
logic, one may conceive it as a proposition about the physical system; a formal statement which 
can be either true or false. 

According to the standard formulation of quantum theory, such an observable is represented by 
a self-adjoint operator A with two distinct eigenvalues a ^ b. Its spectral decomposition has 
the form: 


A = aP + bQ 


(2.41) 
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2 States and effects 


where P and Q are two orthogonal projectors such that P + Q — 1, i.e., P = t — Q. The 
particular values of a and b do not matter as long as they are distinct. In particular, we could 
have chosen a = 1 and b = 0, so that we simply have A = P. We can then interpret that, in 
the outcome of a measurement of P, obtaining the eigenvalue 1 signifies that the proposition 
P is true , and obtaining the eigenvalue 0 signifies that it is false. The converse proposition, its 
negation, is then simply characterised by the observable P L = 1 — P. A state p of the system 
assign probabilities to the truth values of these propositions: 

prob(P is true) = Tr(pP) (2.42) 

prob(P is false) = Tr(pP J ~) = Tr(p(l — P)) = 1 — Tr(pP) (2.43) 


Like in the previous chapter, where we introduced uncertainty about our knowledge about 
the state, one can also add uncertainty to our knowledge of the measurement device. One 
could imagine that, with probability p our measurement device shows the wrong measurement 
outcome. This means that we must review our probability (2.42) as follows: 


prob(P is true) = (1 — p) Tr(pP) + pTi(pP ± ) (2.44) 

= Tr(p((l-p)P + pP ± )) (2.45) 

= Tr(pP), (2.46) 

where we have defined the operator 

E : = (1 -p)P + pP ± (2.47) 

Moreover, (2.43) changes accordingly to 

prob(P is false) = pTr(pP) + (1 — p) Tr(p(l — P)) (2.48) 

= Tr(p(l — E)). (2.49) 


We see that the operator E is not a projector in general, but it is still positive E > 0, and Tr(pP) 
directly yields the probability of the outcome of a measurement. We can therefore generalize the 
concept of proposition , allowing them to be represented by any operators satisfying 0 < E < 1, 
rather than simply projectors. 

These operators are traditionally called effects. It is important to observe that taking the 
spectral decomposition of E would have no physical meaning unlike for an observable. 

Also we note that projectors are a special type of effect, which we call sharp, because they do 
not introduce extra uncertainty beyond that represented by the state. 


2.5 Abstract State/Effect Formalism 

Before we discuss the properties of quantum effects in more details, we introduce the abstract 
state/effect formalism which will set the language for a framework which unifies classical and 
quantum theories. It is known as the state/effect formalism, or generalised probability theory. 





2.5 Abstract State/Effect Formalism 
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Such a theory consists of a set of effects (propositions about the systems which can be true or 
false) and a set of states (an assignment of a probability to each effect). This is essentially a 
Bayesian point of view on physics, where a state represents an observer’s state of knowledge 
about a system. It may or may not represent an objective attribute of a physical system, 
depending on whether the observer is actually right, or whether such objective attributes even 
exist. 

The most basic mathematical structure that one usually requires these sets of effects and 
observables to possess is one that gives the ability to take convex combinations of objects. E.g., 
if pi and P 2 are two states, then one can form a new valid state as pp\ + (1 —p)p 2 where p G [0,1]. 
This corresponds to the physical interpretations of a mixture: the observer is uncertain whether 
the state is p\ or p 2 , and attributes probability p to it being pi. The same can be done with 
effects as in the previous section. 

The following definition characterizes special points of the convex set of states: 


Definition 2.8. 

Let be the set of states. Then the element p G 0 is extremal, or pure, if does not admit a 
proper convex combination, i.e., if there exists 0 < p < 1 such that 

P = PPi + 0--p)f>2, (2-50) 

then one must have Pi = p 2 = p. 


As seen in the previous section, quantum theory in finite dimension can be formulated in this 
language as follows: 


2.5.1 Example: Quantum Theory 


Given a finite-dimensional Hilbert space , we identifies states and effects with certain operators 
as follows: 


states p : p > 0 , Tr(p) = 1 
effects E : 0 < E < 1, 

where a state p assigns to the effect E the probability: p(E is true) = Tr (pE). 

This defines the complete structure of finite-dimensional quantum theory, which is essentially 
the Born rule. We will see below how the other standard postulates, such as “state collapse” 
can be derived from this structure. 

One can check that the extreme points of the set of states are of the form p = if G 

and the extreme points of the set of effects are the projectors: E 2 = E. 

Another example that fits into this framework is classical probability theory: 
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2.5.2 Example: Classical Probability Theory 


Let f] be a finite set. We define states and effects as 

states : p = {pi} i£ n , 0 < p* < 1 , J^Pi = 1 

i 

effects : E = {e;} ie n , 0 < e* < 1, 

where a state p assigns to the effect E the probability: p (E is true) = YliPi^i- 

The states are characterized by a probability distribution on the set If. Extremal states are 
given by p t = S t j , and are hence one-to-one with elements of fh 

Normally, propositions are characterised as subsets co C hi: namely those pure states for which 
the proposition is true. These corresponds to the effects associating 1 with elements of u and 
0 with elements not in u. They are the extremal points of the full set of effects. 


For example let If = Z. If this system characterises a “random variable” called X, then the 
proposition U X is positive” corresponds to the effect E = {eff i&L with 


e* = 



if i > 0, 
otherwise, 


(2.51) 


or, more generally, the proposition ‘W belong to some subset u C hi” corresponds to the effect 


e* = 



if i E u , 
otherwise. 


(2.52) 


Evaluating the probability for the latter amounts to 

P (E is true) = p(X E u) = e iPi = ^2 Pi, 

i i£uj 


as expected from standard probability theory. 

This framework allows us to formulate classical probability theory in terms of a quantum theory 
with additional restrictions: Let {pffi^n denote a classical state and E = {ej}j e Q a classical 
effect on an outcome set fl Enforcing that the states p and Effects E are diagonal in a fixed 
basis {|i )}” =1 allows for the following identification: 

P = Y,Pi\ i )( i \ 

i 

E = l*X*l- 

i 

Therefore, this amounts to valid state/effect configuration of a quantum system, where the 
probability assigned to the effect is exactly the same as in the classical setting p (E is true) = 
Tr (pE) = Y.iP%e%. 
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2.6 Algebraic Formulation 


The state/effect formalism allows us to introduce a formalism which generalises both quantum 
and classical theories, while keeping as much structure as possible from both. Moreover, this 
framework allows for infinite-dimensional generalizations even beyond Hilbert spaces, though 
in this lecture we focus mostly on finite systems for pedagogical reasons. 

For this purpose, recall the notion of an algebra A over C. ft is a complex vector space on which 
a multiplication is defined. This multiplication is compatible with addition, i.e. a(b+c ) = ab+ac 
for a, b, c G A. The multiplication is called associative if also a(bc ) = a{bc). Also, it is called 
unital if it contains a unit 1 for the multiplication: la = al = a. Based on these structures we 
define 


Definition 2.9. 

A ^-algebra A is an associative algebra with an anti-linear map * : A —> A satisfying the 
following properties Va, b G A: 

. (a*)* = a 

. ( a &)* = b*a* 

Moreover it is unital if it contains a unit 1 such that, Va G A, la = al = a. 


A prominent example of a ^-algebra is given by the complex matrices where the involution * 
is given by the conjugate-transpose operation f. Of course, there is a natural concept of 
substructure of ^-algebras: 


Definition 2.10. 

Let B be a ^-algebra and A <Z B. A is a *-subalgebra, if the following properties hold: 

• A is a linear subspace: a, b G A, a, /3 G C =>■ aa + f3b G A 

• A is closed under multiplication: a, b G A =>■ ab G A 

• A is closed under conjugation: a G A =>■ a* G A 


The structure of a ^-algebra allows us to introduce a concept of positivity that generalises that 
of operators: 


Definition 2.11. 

Let A be a ^-algebra. We say a G A is positive, i.e. a > 0, if there are elements bi G A , 
i — 1,..., n such that 

a = J2 b i b i- 

i 
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This definition is such that the set of positive elements is a cone, namely closed under linear 
combinations with positive scalars. In turn, this makes the relation a < b defined by b — a > 0 
into a partial order. 


It is not hard to see that for matrices, with A* := A', this condition is equivalent to Dehni- 
For instance, if A is positive as a matrix, we have seen that A = Y^i a i |*)(*|, where 


tion |2.2 
a* > 0. Hence, using Bi := 


\ 'l 


we have A = Yhi Bj B,. 


Using this structure, we can associate a generalized probability theory to any unital ^-algebra 

A: 


Definition 2.12. 

Let .A be a (possibly infinite-dimensional) ^-algebra. We define the effects associated with A 
to be the positive elements which are smaller or equal to the identity 

effects: e E A, 0 < e < 1 

Moreover, we define the states as certain linear functionals on A: 

states: / : A —> C, / linear, positive and unital, 

where positivity of / means that /(a) > 0 for all a > 0, and unitality is the normalisation 
condition /(l) = 1. 


Note, that this assignment gives a physical meaning to the abstract mathematical object of a 
*-algebra, once some effects are associated with actual experimental setups. 

Before we describe examples, we need the following correspondence between linear functionals 
on some finite dimensional algebra and elements inside the algebra. 


Theorem 2.2. Suppose that A is a *-subalgebra of the n-by-n complex matrices Ai n : A C M n . 
Then for every linear functional f : A —» C there exists R e M n ( C) such that 

f{A) = Tt(RA) \/AeA (2.53) 


Proof. Since 1 = |i)(i|, we have A = Eij{i\A\j)\i){j\. Therefore, f(A) = Eij(i\A\j) f(\i){j\). 

If we define R := 'Ey f{\i)(j\)\j)(i\, then f(\i)(j\) = (,j\R\i), and f(A) = Jfij{i\A\j)(j\R\i) = 
Tr(RA). □ 

Hence, the dual (that is the the set of linear functionals) of a finite dimensional algebra is 
isomorphic (as a linear space) to the algebra. This fact does not translate into the infinite 
dimensional setting. 


The previously defined probability theories can in fact be defined in this way from ^-algebras: 
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2.6.1 Example: Quantum Theory 


Given a Hilbert space H = C d , consider the unital ^-algebra of linear maps on R: A = 13(H) 
(B stands for “bounded linear map”, which is true of any map in finite-dimension). In this 
case, the star operation is given by the adjoint: 

A* = A*. 

This algebra is the same as that of d-by-d complex matrices: A = Ad,/. Building the states and 
effects from A yields quantum theory with Hilbert space H. 

Indeed, the set of effects E is precisely the one we previously identified in quantum theory. 
Therefore, we just need to show that the states defined as positive linear functionals / : B(T-L) —> 
C correspond to density matrices. Theorem |2.2| tells us that for every linear functional / we 
can find a Re Aid such that f(A ) = Tr(/L4) \/A e B(H). The positivity of the functional / 
says that for all positive matrices A > 0, Tr(RA) > 0. Choosing in particular A = |'0)(-0|, this 
says that ( , 0|-R|V’) > 0 for all ^ e H, i.e., R is a positive matrix. The trace-condition Tr R — 1 
also follows directly from /(l) = 1. Therefore R is a density matrix. 

Likewise, we also obtain classical probability theory if we consider a commutative ^-algebra: 


2.6.2 Example: Classical Probability Theory 


Given a finite set 12, we consider the unital ^-algebra 

A = L(12) 

of functions a : 12 —> C, with all algebra operations defined pairwise: for any a, b G A, x G 12, 
a/3 G C, 


(era + f3b)(x) := aa(x) + /3b(x ) (2.54) 

(ab)(x) a(x)b(x) (2.55) 

a*(a;):=a(x) (2.56) 

where a denotes the complex conjugate of a G C. 

It is easy to check that positive elements of A are those for which a(x) > 0 for all x , and 
hence, the effects are precisely those of classical probability theory. A natural linear basis 
of A is given by the index functions Xx defined by Xx(y) — $x, y f° r ah x,y € 12. Hence, 
given a linear functional / : A — » C, we have /(a) = Y^ x a i x )f(Xx)- Let p e A defined by 
p(x) = f(Xx), we have /(a) = J2 X a (x)p(x). If / is a state, then we must p(x) > 0 for all x and 
/(l) = 12 x p(x) = 1- Hence states are probability distributions over 12. 

Observe that this algebra is equivalent to that of diagonal d-by-d matrices, where d is the number 
of elements in 12. Indeed, if we identify any a € A with the diagonal matrix A with diagonal 
elements a(x), x G 12, which we can also write using the bracket notation A = Ylx^n a(^)|^)(x|, 
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where the vectors x), form an orthogonal basis, we see that the matrix product and other 
operations are the ones above. 

We will also need to consider classical systems with continuously many pure states. This can 
be formalised in several ways. Here, we will represent this situation by a measurable set fi, i.e., 
a set with a concept of volume, or measure, for most subsets. A typical example will be = M. 
The algebra that we will consider, 

L(fi) = L°°(fi), 

is the commutative algebra of bounded functions / : O — > C, i.e., there is a constant C such 
that |/(x)| < C for all x G except possibly for a subset of measure zero. L°°(fl) is the natural 
classical counterpart of B{H). 

In addition to the ^-algebra structure described above, one normally also requires a norm || ■ || 
to be defined (and finite) on every element of the algebra, and which respects compatibility 
conditions with respect to the product and * operations. 

If also ||x?/|| < ||x|||M| for all x,y e A, then A is a Banach algebra. If, moreover, ||x*x|| = 
||x|| ||x* || for all x £ A, then A is a (W-algebra, which is the most popular algebraic framework 
for quantum theory. We define here this important concept for the reader’s benefit, but we will 
continue working with concrete operators on Hilbert spaces for simplicity in this document. 


2.6.3 Example: C*-Algebra 

A C*-algebra is a ^-algebra that is also a Banach algebra, i.e., complete with respect to a norm 
satisfying \\xy\\ < ||x||||?/|| for all x,y £ A, and also such that 

11xx*|| = ||x||||x*|| = ||x|| 2 for all x E A. (2-57) 


A typical example is given by the set of bounded linear operators BifH) on a Hilbert space BL, 
where the norm is the usual operator norm given by the square root of 


Pir = sup 

if>en 


(iJ)\A* A\iJ)) 

(?/#) 


We see that the requirement of a finite norm eliminates unbounded operators such as most 
quantum mechanical observables (position, momentum, energy). However this is not a problem 
, since, as seen above, the interpretational content of these observables can be replaced by 
elementary effects (in this case projectors) which are always bounded. We will come back to 
this point later when we need to talk about infinite-dimensional Hilbert spaces. 


The above axioms can be justified by the fact that they precisely capture the nature of algebras 
of linear operators: 


Theorem 2.3. Any C*-Algebra is a norm-closed *-subalgebra of BifH) for some Hilbert space 

H. 
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The proof is beyond the scope of this document. Instead, we will focus on concrete ^-algebras 
of matrices. 


2.7 Structure of matrix algebras 


In the following we would like to characterize the general structure of the “matrix algebras”, 
namely the *-subalgebras of the algebra of complex matrices of a given dimension. Because of 
the finite-dimensionality of these spaces, these subalgebras are also automatically C*-algebras 
with the operator norm ||H|| defined by 



max 

b 


(if\A^A\if) 

(V#> 


However, contrary to the infinite-dimensional case, we will not need to make use of the norm 
for characterising these algebras. 

The smallest such algebra that one can build is the algebra alg(H) generated by a single self- 
adjoint matrix A, i.e., all linear combinations of (non-zero) natural powers of A. This algebra 
is commutative and unital, and can be characterised as follows: 


Theorem 2.4. If A is a self-adjoint matrix with spectral decomposition A = E,; a l P l (where 
at jtz 0 are distinct eigenvalues of A, and Pi, - ■ ■ , P m are orthogonal projectors), then 

alg(H) = span (Pi,..., P m ) 


namely, the *-algebra generated by A is equal to the span of the spectral projectors of A. 
plicitly, 


Pi = n 

*¥=3 


A - ail A 

aj — ai 


Ex- 


where 1 A = A 0 = E™ i Pi is a unit for the algebra. 


Proof. We used the notation 


alg(H) 



and span(Pi,..., P m ) 



It is clear that any power of A is inside this space spanned by the projectors Pi, since A n = 
Ej afPi- What we have to show, however, is that the projectors P* can be written as polynomials 
in A. But first, we need to show that alg(kL) contains the unit element A 0 . Since the set of 
eigenvectors of A is complete, we have Db( a il — A)A = 0. Indeed, this polynomial sends every 
eigenvector to zero. Multiplying by the generalised inverse of A we obtain bb; ( ai A° - A) = 0. 
Since a* 7 ^ 0 for all i, we can simply solve for A 0 . We can then obtain the projectors explicitly 
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as linear combinations in A 1 for i — 0 ,..., n: 

(cifc <2.1) (Ofc dj—l) (®fc Qj+l) ®n) j-) 

tj = 2 ^-• • •-• • •-“fc 

CL d £Zji (Zj — d (Zji CL j _|_ d ZZ^ 


°k,j 


E 


(flfci - «l) (%-i - a?-l) (flfcj+1 - %+l) (flfcn - On) 


A*1 A* • k ■ . 1 A- ^ 7 ^1 


(Z j d zz ji CLj _|_ d 


CL j ZZ^ 


P P P 

1 ki ■ ■ ■ 1 kj-i 1 k j+ 1 • • 


^k 1 ,k 2 ^k 2l k 3 r-- 


E 


(G'/ci ^l)-P&i iflkj — i Q'j—l)Pkj—i {flkj -|-i ^' < 7+l)-Pfc J _|_i (^fc n Q*n)Ph 


A*1 k 1 k < 1 k O'! ^1 


a j a j -1 


a j a j +1 


dj CL n 


TT ® i)Pk _ tt X)fc ®fc-Pfc Pk _ tt ^4 

~r a,- — a, ~r. " ■*■■*■ 

* 7=7 j 1 i/j 


ZZ 


a j a * 


□ 


This result has a natural generalisation to arbitrary matrix algebras: 

Theorem 2.5. Hn?/ *-algebra A of complex matrices is unital, and is equal to the span of the 
projectors it contains: 

A = span{P G .4. : P^P = P}. (2.58) 


Proof. If we consider an arbitrary element B £ Awe can write it as combination of a hermitian 
and anti-hermitian part, i.e. 

B + B* B - B* 

B =--- H -. 



=:Im (B) 

Re (B) and Im (B) are both self-adjoint by definition, and manifestly inside of the same *- 
algebra A. Moreover, theorem 2.4 showed that the spectral projectors of these matrices are 
also within A, since it must contain alg(Re(B)) and alg(Im(B)). Hence if we write their 
spectral decompositions as Re (B) = Xa biPi and Im (B) = Xa biQi then we obtain 


B = E a i p i + 1 E b jQj, 


where Pj, CP G A. This proves Equ. 2.58 In order to build the unit element, let us write simply 


Ti, i = 1,..., n for the finite set of projectors including the Pf s and Qf s that we just built. 


These matrices span the whole of A as we have just seen. From theorem |2.4[ we know that the 
projector P := (XaP° on the range of X)?: T % is also inside A , and satisfies P(X,; Tf) = Jf , T % , 
which implies Xa(l — P)Tj(l — P) = 0. But this is a sum of positive operators, therefore 
each one of them must be zero: (1 — P)7)(l — P) = 0 for all i, as well as their square roots 
Tj(l — P) = (1 — P)T t = 0 (which can be proven by taking the expectation value of the previous 
expression in an aribtrary vector). It follows that T*P = PX) = Ti for all i. Since the matrices 
Ti span A, this shows that P G A is the unit of A. □ 


Note that Equ. 2.58 is also true in any ^-algebras of operators on a Hilbert space, even an 
infinite-dimensional one, provided that the algebra is closed in the weak operator topology , i.e., 
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that if a sequence of operators A n are such that {ij}\A n \ijf) converges to ('if\A\'i() for all vectors 
\if) then A is also inside the algebra. These algebras are special types of C*-algebras called von- 
Neumann algebras. The span of the projectors must then also be closed in that topology. There 
are, however, many non-trivial C*-algebras which contain no projector besides the identity and 
zero elements. 

The unital algebra alg(kL) generated by a single self-adjoint matrix A is the prototype of a 
commuting matrix algebra. Indeed, all commuting algebras are of this form: 


Theorem 2.6. Any commuting *-algebra A of complex matrices is of the form 

A = alg(H) = span (Pi,..., P n ), 

where A is a self-adjoint matrix and Pi ,..., P n , are a complete set of orthogonal projectors, 
and n is the dimension of A (as a linear space). 


Proof. Consider a linear basis Bi,..., B n of A. By decomposing each Bj into a real and 
an imaginary part, one obtains a family A\,... ,A m , where m = 2n, of self-adjoint operators 
spanning A. Now consider the spectral decomposition of each of these operators: Aj = , ajQj. 

Since the algebra is commutative, these projectors also commute with each other. This implies 
in particular that the product of any two such projectors is also a projector. Indeed, if P 
and Q are two commuting projectors, then (PQY(PQ) = QPPQ = QPQ = PQQ = PQ. 
Consider therefore the projectors Pi = Q\ x ■ ■ ■ Qf , where G {1,..., m}. These projectors 
are all orthognal to each other. Indeed, if / and J differ by a single entry, then Pi and Pj 
contains in their product two orthogonal projectors, which, given that they all commute, implies 
PiPj = 0. Moreover, the span of these operators include also the projectors Pf , and hence span 
the algebra. Hence we have obtained a complete family of orthogonal projectors Pi spanning 
A. Observe, however, that most of them are typically equal to zero. Indeed, due to their 
orthogonality, they are linearly independent, and hence only n of them can be non-zero. Let 
Pi,... ,P n be those nonzero projectors. If we then pick a distinct real number a t for each i, 
then the matrix A = a iPi generates the whole of A (since each Pi is inside alg(H) Cd). □ 

Now that we have understood the form of commutative matrix algebras, we gain a first handle 
on generical matrix algebras by considering a commutative algebra that they contain; their 
center: 


Definition 2.13. 

Let A be a ^-algebra. The center of A is the commutative algebra 

Z{A) = {A e A : [A, B] = 0 , VB e A}. (2.59) 


Hence the center Z(A) is made of those elements of A which commute with all the other 
elements of A. Clearly, Z(A) is a ^-algebra: If we take A, B G Z(A) we see that also AB G 
Z(A). Indeed, for all C G A, ABC = ACB = CAB. Moreover, Z(A) is clearly commutative 
(abelian), i.e., for all A, B in Z(A), [A, B] = 0. 
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The structure of the center Z(A) gives us already important information about the structure 
of a general matrix algebra A. In order to express the resulting property, however, we need to 
introduce the concept of direct sum © for algebras. The algebraic direct sum © can be defined 
concretely as follows: given two square matrices A and B , respectively of sizes nx n and m x m, 
we define the matrix A © B of size nrn x nm as 


A®B 




where the 0’s represent rectangular matrices of the apropriate sizes with only zero components. 
It should be clear how this generalises for the direct sum of more than two matrices. 


This direct sum could also be defined algebraically also on general abstract algebras A and B: 
A © B is the direct sum of A and B as linear spaces equipped with the product defined by 

(A © B){A! © B') := (AA') © ( BB'). 

As block matrices, this is simply 

(A 0\ (A' 0 \ _ /AA' 0 \ 

\0 BJ\ 0 B'j ~ ^ 0 BB') ‘ 

With the help of this concept, we obtain the first importance piece of information about general 
matrix algebras: 


Lemma 2.1. For any *-algebra A of complex matrices, there is a unitary matrix U , and smaller 
matrix algerbas Ai with Z{Af) = Cl, such that 

UAU ] = Ax © ■ ■ ■ © A n © 0. 


Proof. Theorem 2.6 tells us that there exists a minimal set of orthogonal projectors {Pi}f =1 
such that Z(A) = span (P 1 ,..., P n ). Note that the sum of all projectors is the identity of A, i.e. 

Pi = 1 a- This matrix 1 A should not be confused with the identity matrix: it ould be in 
general any arbitrary projector. 


Due to the fact that the Pi commute with every element of the algebra, we obtain the following 
decomposition of any element A G A: 

A = 1 a A = ( ^ Pi ) A = PiKA = Y, p i Ap i ~ A i © A 2 © • • ■ © A n © 0, (2.60) 

for some smaller matrices A* whose dimensions sum to the total dimension of the projector 1^. 
The symbol ~ means that the two matrices are equal up to a change of orthonormal basis. 


In this basis, we have, 


1a — Idi © ■ ■ ■ © 1 d n © 0 rfn+1 
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where 1^ denotes the identity matrix of dimension di, where d, t = Tr P t , and 0 d n+1 is the zero 
square matrix of dimension 

n 

dn+ i • d ^ ) di 

i =1 

where d is the dimension of the original matrices in A, such as A. This is just the identity 
matrix accept for a few trailing zeros on the diagonal ( d n+ \ of them). 

Also, the smaller matrices Ai are defined as representing the non-trivial parts of the matrices 
PiAPi in the new basis: 


P,AP, ~ 0 




Ai © 0 d. 


'i+iH- \~d n +1 * 


Hence, in this special basis we have, in block matrix notation, 


( Al 

tl 2 


A ~ 


V 


\ 



Ai © A 2 © • • • © A n © 0. 


where the blank spots are filled with zeroes, and 0 represents the zero square matrix of dimension 
dn +1 • 


Because the projectors P* only depend on the algebra A, not on the chosen element A, all 
elements of A have this particular form when expressed on that specific basis (but for different 
values of the small matrices Ai). In fact, it is easy to see that the allowed values of the matrices 
Ai must also form a ^-algebra that we call A t , so that we obtain: 

A ~ Ai © • • • © An © 0. 


□ 

In order to completely elucidate the structure of A, we still need to understand the structures 
of the smaller algebras A%. What distinguishes them from a general matrix algebra like A is the 
fact that their center is “trivial”, i.e., it consists of all multiples of the identity: ZfAf) = Cl. 
Indeed, if there were other elements in the center of Ai, there would be a hirer decomposition 
of the minimal set of projectors of the center. Such a matrix algebra with a trivial center is 
called a factor. 

We now show that such a matrix algebra factor always consists of all matrices of the form 1 © A 
in some basis: 

Lemma 2.2. Let A be a *-algebra of d-dimensional complex matrices, with trivial center (i.e., 
a factor): Z(A) = Cl^. Then, there is a unitary matrix U such that 

UAU ] = t m ®B(C n ). (2.61) 


with mn = d. 
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Proof. Our proof is adapted from Ref. [5j. Given any nonzero element A G A, consider the set 
X A := {XAY : X,Y G *4.}. The set Z A is an ideal , i.e., for all B G A and A G Z A , AB G Z A 
and BA G Z A . Let us show that Z A contains the identity matrix 1 (which implies Z A = A). 
It is clear that Z A is closed under multiplication and under the * operation, and hence forms 
a matrix algebra. From theorem 2.58 we know that Z A contains a unit P G Z A . This implies 
that for all B G A, P{BP) = {BP) = {BP)P since BP G Z A . Since also (BP)P = BP (P 
is a projector), we have PB = BP and hence P G Z{A). But since A is a factor this implies 
P = 1 G Z A . 


In other words, for any A G A, 1 = Yli.XiAYi for some X tl Y i: G A. Hence, for all A, B G A, 
B = JfiBXiAYi. If B 0, this implies that at least one of the terms in the sum must 
be nonzero: BXiAY t ^ 0, which implies BXiA ^ 0. We have therefore shown that for all 
A, B ^ 0 G A there exists X G A such that BXA ^ 0. Although this may not seem very 
impressive, this is the essential property of A we need. 

Now let us consider a maximal commutative *-subalgebra C of A. The maximality means that 
if B G A is such that [B, A] = 0 for all A G C, then B G C. This exists because we can build it 
by progressively adding any such element and its dual B^ to a given commutative subalgebra 
until there is not more left. 


From theorem 2.6 we know that C = span(Pi,..., P n ), where Pi are a complete family of 


orthogonal projectors. We have just shown above that for every pair i,j, there exists a X VJ G A 
such that 

Fij '■= PiXijPj ~f~ 0 . 


Observe that F} 3 F t] commutes with every elements of C and is therefore contained in C due to 
its maximality. This means that it is a linear combination of the projectors Pf,, but its product 
with Pfc, k ^ j is zero, hence 


/' /' . — P X PX P = a- P 


for some a t j > 0. In particular, this implies that all the Pf s project on spaces of same dimensions 
m since F- 3 F t j above maps all vectors in PjC d via PiC d back to PjC d without sending any to 
zero. Therefore, we can decompose C d into C m 0 C n with an orthonormal basis | i) 0 | j) such 
that 

m 

p i = J2 \j )01 ® l*)(*l = 1 81 |i)(i|. 

3 = 1 

This yields F\j = Y)j0|i)(j|, where Yij = (l0(i|)Wj(l0|j))- But then F^Fij = Y,^Yij^\i){i\ = 
a tJ t 0 \i)(i\, which implies Y^Y l3 = a^t, and hence 

F^ = Oij 1 0 \i)(j\. 


We only have left to show that the matrices F t j span A. Consider any A G A, for the same 
reason as above, we have PiAPjXjiPi oc. Pi. It follows that 

PiAPj = PiAPjPj oc PiAPjF^Fji oc oc P % X l3 P 3 . 

This implies that there is are linear functionals fij : A —>• C such that, for all A G A, 


PiAPj = fij(A)F tj . 
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Therefore the matrices span A, since YiPi — 1 and hence A = Yij PiAPj = Yij fij(A)Fij. 


O' 


It is now straightforward to combine Lemma |2.1| and |2.2| in order to completely elucidate the 
structure of general matrix algebras: 


Theorem 2.7. For any *-algebra subalgebra A of B(C d ), there is a unitary matrix U such that 


UAU ] 


r n 


0 t mi ® M ni (C) 

. 2=1 


© Orf 0 


(2.62) 


where YiLi n i m i + d 0 = d. 

Equivalently, this means that the elements of A are precisely those matrices of the form 

(t mi ®Ai \ 


A = U ] 


Imjv © A n 


V 


u 


0 / 


for any set of square matrices Ai of size ni. 

In this notation, the center Z(A) of A is made of all matrices of the form 


( OL\ Imi © Ini 


A = U ] 


© Injv 


V 


\ 


0 


u 


for any family of complex numbers o-i,... Qn- 
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(2.63) 


(2.64) 


Theorem 2.7 allows us to understand the type of physical system represented by a generic 
matrix algebra A, in terms of the state/effect formalism. For short, we write 


A — 1 ® Ai 


(2.65) 


for generic element of A, according to the decomposition given by Theorem 2.7 For conciseness, 
in this notation we leave implicit the possible trailing zero block. Also, one must remember 
that this is block-diagonal only in a basis which may not be the canonical one. 


It is easy to see that A^ — A if and only if A\ = A, for each i. Moreover, by writing A in 
diagonal form, we immediatly see that the eigenvalues of A are all in the interval [0,1] if and 
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only if that is the case also for each A- Therefore, we conclude that the effects 0 < A < 1 of 
A are precisely the operators of the form © ; 1 0 Ai where 0 < Ai < 1 for each i. 


From theorem 2.2, we know that the states can be represented by matrices P as the functionals 


A m- Tr(PA) 


on effects A. But since the set of effects is restricted, two matrices R and R' may actually 
represent the same functional, i.e., Tr(iL4) = Tr(P'A) for all A e A. 


Let us find a special matrix R giving a unique representation of the functional. Let Pi,..., P/v be 
the projectors spanning the center of A. We know that the elements of A satisfy A = J2i RAP t . 
Therefore, 

Tr (RA) = E TRRRAPi) = '£ i Tr(P i R P i A). 

i i 

This means that the functional is uniquely specified by the matrix 

R! = E p i Rp i = ® R > - 

i i 

There is however still some ambiguity left on the matrices P* since 

E Tr(Pi(l 0 A)) = E Tr ((! « Tr'(P i ))(l 0 A)) • 

i i 

We obtain that the same functional is uniquely represented by a matrix of the form 

R = (J) 1 0 Ri, 


namely an element of the algebra A itself. 

Using the same argument, as for the effects, this matrix R is positive if and only if R { > 0 for all 
i. The last condition for R to represent a state is the normalisation condition TrP = 1. Since 
Ri is positive, we can write Ri = riPi where pi is a normalised density matrix and r t := TrP*. 
The normalisation condition then says 


Tr P = E m i r i — 1- 

i 

Let us write p, = miTy. These numbers form a probability distribution, and we have 

R = ® Pi — ® Pi, (2-66) 

i m i 

which is a unique representation for a state of the theory defined by the algebra A. 

We make two observations. The first is that P is also a density matrix on the quantum system 
defined by the full matrix algebra of which A is a subalgebra. Therefore, restricting the quantum 
effects to those of the subalgebra A is equivalent to restricting the density matrices also to A. 

The second observation, is that a state of the system defined by A is represented by a proba¬ 
bility distribution pi, ... ,Pn, hence a classical system , together with a set of quantum states 
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Px ,..., pn • Hence it is a sort of hybrid quantum/classical system. This interpretation is justified 
by the fact that these states characterise knowledge about the outcome of a procedure where, 
say, a different quantum state is prepared depending on the value of a classical variable. 


Moreover, it is clear that if A7 = 1 , then M is just a quantum system, while if the dimension of 
each factor Ai is equal to 1 , then this is just a classical system. 


A particular example of such a hybrid system is one which consists in a classical system next 
to a quantum one. Within the algebraic framework, we can compose systems of different type. 
For this purpose, we generalise the tensor product introduced in Section 2.3.2 in the context 
of quantum systems to systems defined by arbitrary algebras, Mi, M 2 —» Mi (8) M 2 . Algebras 
are vector spaces, for which the tensor product is defined in the usual way. In addition, we 
must defined how the product and * operations, which are specific to the algebraic structure, 
are built on the composed systems. Specifically, we demand the following properties for the 
multiplicative structure: 


(Mi 0 M 2 )(f?i 0 B 2 ) — A\B\ 0 A 2 B 2 (2.67) 

(Mi 0 M 2 )* = A{ 0 A* ( 2 . 68 ) 

Importantly, there are canonical embeddings of a single algebra into the tensored space, which 
allow one to think of each algerba A % as a subalgebra of Mi 0 M 2 , via the homomorphisms 


Mi —y Mi 0 M2 M2 — y Mi 0 M2 

M M 0 1 M 4 I 0 M 


In case of matrix algebras, which are *-subalgebras of full matrix algebras, the tensor product 
is precisely the same as the quantum one, which is also called the Kronecker product, e.g., 


'Mu • 



^ dwB 

CLi n B\ 



®B = 




^nn / 


\CLn\B 

UnriB J 


where a l3 are the components of a matrix, and B represents another matrix. 

For instance, if we compose two finite dimensional quantum systems Mi = 13(Hi) and M 2 = 
B(T-L 2 ), the tensor product is isomorphic to bounded operators on the tensor products of the 
hilbert spaces Mi 0 M 2 = B{hi\ 0 H-i). 

Recall that a classical system associated with the finite set of pure states (sample space) 
corresponds to the commutative algebra of functions M = L(Q) on Q, namely the space of 
vectors with elements indexed by elements of fi, equipped with the component-wise product. 
Taking the tensor product of two classical systems Mi = L(Q 1 ) and M 2 = L(Vt 2 ) results in 
the set of functions on the cartesian products of the sample spaces: Mi 0 M 2 = L(Q 1 x f^)- 
Alternatively, one may also view those commutative algebras as sets of diagonal matrices, with 
diagonal elements index by H. The Kronecker product then gives the same result. 

A natural hybrid system then is obtained by considering the composition of a classical system 
Mi = L(Q) and quantum system M 2 = B(H). The easiest way to see what happens is, again, 
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to represent L(Q) as a set of diagonal matrices. Using the Kronecker product, one then find 
that the tensor product yields a block-diagonal algebra 

where each factor B(7~L) on the right hand side of the equation is indexed by an element of fh 


For instance, a completely uncorrelated state on A\ 0 A 2 has the form 



r p\ 

\ 


'PiP 

\ 

III 

0 

y 



0 p = 




\ 

Pn) 


\ 

PnP) 




where the orthogonal vector states |i) just serves for the matrix representation of the classical 
state {pi}fLi = J2iPi\i)(i |- 


Similarly, a general correlated classical/quantum state has the form 


(PiPi 


V 


\ 


PnPnJ 


N 

1=1 


pi 


This is a special case of the general expression Equ. (2.66), where |f)(z| is simply the identity 
operator on a one-dimensional factor, and where each state p, leaves in a Hilbert space of same 
dimension. 
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We want to identify the most general way of representing a transfer of information from a 
system represented by some algebra A\ to another system represented by A^- This can be 
done in two ways: either one maps states to states (Schrodinger picture), or effects to effects 
(Heisenberg picture). I 11 infinite dimensions, the Heisenberg picture is more general. However, 
we begin with the Schrodinger picture because it is somewhat more intuitive. 

Mathematically, a map preserving all the structure of A\ is given by the following definition: 


Definition 3.1 (^-homomorphism). 

Let A\ and A 2 be ^-algebras. A map 0 : A\ —> A 2 is referred to as ^-homomorphism (or 
^-algebra morphism) if it satisfies 

a) (j) is linear 

b) <j>{AB) = <t>{A)<t>{B) 

c) (j>(A*) = 0(A)* 
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However, these conditions are too strong for our purpose. Indeed, we only want to preserve 
the state/effect structure. This structure does rely on the algebra product, but only through 
the partial order defined via the concept of positivity. Indeed, recall that positivity is defined 
using the product (Definition 2.11). But, as we will see, the requirement that a map preserves 
this partial order is much weaker than requiring that it preserves the product structure. 


3.1 Schrodinger Picture 


A map £ on density matrices should preserves convex combinations of states. Indeed, since 
the state J2iPiPi —where {p;}f =1 is a probability distribution and pi are density matrices— 
represents the belief that the state is p t with probability p l , its image under the map £ should 
be 


namely, the image of p under £ ought to be state £(pi) with probability /y. Since the set 
of density matrices for the system represented by A\ span the whole of Ai, this map can be 
naturally extended to a linear map 

£ \ A\ —y Ao 


Moreover, for the image of any density matrix to also be a density matrix, the map £ must 
send positive matrices to positive matrices (according to which we simply say that £ is itself 
positive ), and it must also preserve the trace of matrices: 

X > 0 =y £(X) > 0, 

and 

Tr(£(X)) = Tr(AT). 

But there is an extra more subtle consequence of the positivity condition. Namely, £ should 
also respect all those conditions if we see as part of a larger system. That is, if we add any 
system B{T~C) on which £ acts trivially. On the large system, this is represented by the map 
(£ ® id) : A\ Cg) B(%) —y A 2 0 B(di) defined by 

{£ <g> id){X ® Y) = £{X) ® Y. 

Classically, one would expect that the positivity of £ implies that of £ ® id. Indeed, if the 
algebra A\ is commutative, then a general operator in A\ ® B{l~i) is of the form 

i 

where X { 6 B(H). Since the eigenvalues of A" are just those of the A^’s, it is easy to see that 
X > 0 if and only if X % > 0 for all 1 . We then have 

(£®id)(X) = ^^(|*)(*|)0W. 

i 

Since £{\i){i\) and A,- are positive, so is the right hand side of this equation. This shows that, 
when A\ is commutative, and hence represents a classical system, the positivity of £ implies 
that of £ Cg) id. However, the following example shows that this fails when A\ is quantum: 
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3.1.1 Example: Partial Transpose 

Consider the transpose map £(p) = p T , where p is state on 1-La- Since the transposition 
leaves the eigenvalues of p invariant, this is clearly a positive map (and trace preserving). 
However, consider the setting if we compose the Hilbert space with an ancillary space of the 
same dimension d: H — Ha ®Hb- We extend S by the identity on this ancillary space and 
consider the action on the maximally entangled state |fi) = |*)a <E> J i)b' 

(£ ® i)(|H)(H|) = (s ® i) (^ Y, I*)01 ® l*)OI) = \ Y, OX*I ® l*)0l =: p ( 3 - 69 ) 

\ ij ) ij 

But the operator p is not positive. Indeed, consider the case d — 2, and the state |^) = 
(101) — 110)). A direct calculation shows that 

PW) = -W), 


from which it follows that (i(;\p\i[}) = —1. 

We conclude that positivity is not a sufficient criterion for a map to fit into the state/effect 
formalism of quantum systems. We therefore need to extend this notion to that of complete 
positivity: 


Definition 3.2 (Complete Positivity). 

A linear map £ : Ai —> A 2 is completely positive , iff £ (8) id n : Ai (8) B( C n ) —>■ A 2 <S> B(C n ) is 
positive for all n, i.e. 

Vn G N : (£® id n )(X) > 0 VX > 0. (3.70) 


This gives us all the tool to define a general channel in the Schrodinger picture: 


Definition 3.3 (Channel). 

A channel represents a general transfer of information from a system represented by that algebra 
A\ to that represented by A 2 - When those are matrix algebras, it can be represented by an 
arbitrary linear, completely positive, trace-preserving (CPTP) map £ : Ai —> A 2 , meant to be 
applied to states (Schrodinger picture). 


In what follows we give some simple examples of channels. 


3.1.2 Example: Homomorphisms 


A *-homomorphism (j) : A\ A 2 is & channel if it is trace-preserving. In order to show 
this, we just need to show that it is completely positive. Moreover, since d 0 id, is also a 
^-homomorphism, we only need to show that an arbitrary ^-homomorphism like (j) is positive. 
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This can be done by using Definition 2.11, which says that A > 0 if and only if A = E*Ei. 
We then have 

<KA) = !>(£?£,) = Y.'KEiT’KEi). 

i i 

The right-hand side is manifestly positive, which completes the proof that 0 is a channel. 


An example of ^-homomorphism is that induced by a unitary operator U (defined by the 
property WU = UU 1 = 1), through 


m = u p ia 

In fact, we will see below that when A\ = A 2 is a matrix algebras, then a ^-homomorphism is 
necessarily of this form. 


3.1.3 Example: Mixture of unitaries 


The set of channels between two given algebras is convex. Indeed, if £\ and 82 are two channels, 
and p a probability, then we can define the convex combination £ = p£ 1 + (1 — p)£ 2 via 

£{p) = p£ 1 {p) + (1 ~p)S 2 (p). 

This map £ is clearly trace-preserving. It is also easy to check that it is completely positive, 
because £ 0 id is a convex combination of £\ 0 id and £2 0 id , and the convex combination of 
positive operators is still positive. 

This gives us a first example of channels which are not homomorphisms, and hence go beyond 
unitary maps: Let A\ = A 2 = be the input and output algebra. Take Ui : V. —» V. be 

a set of unitary linear operators and {pi}^ =1 . Define the map £{p) = Y^i PiUipU-. It is trace¬ 
preserving and completely positive. This channels describes the situation in which the unitary 
Ui is applied with probability i. Indeed, the resulting state is the convex combination of the 
possible outcomes, which is compatible with the interpretation of the convex combination of 
states. 


3.1.4 Example: Classical channels 

This definition of channel applies to any system defined by an algebra, not only quantum 
systems. Let the input and output algebra of our system be classical Ai = L(Di),^4 2 = L(£l 2 )- 
We have seen above that whenever Ai is classical (commutative) then positivity is sufficient. 
The linearity of £ makes it possible to write the action of £ on a probability distribution 

{Pi}* e Oi clS 

£(fe}) = {Vi = J2 7T ijPj}iett 2 
0 

The positivity of £ implies nij > 0, V«, j. The trace-preserving property amounts to 

E * = E % = E Pi V {Pi )f=i 

i ij j 
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Choosing the particular distribution pj = 5ij we obtain J2iEij — 1 Vj. This means that a 
classical channel is represented by a stochastic matrix n which is the well-known concept of 
channel from usual probability theory. 


3.2 Heisenberg Picture 


Consider the algebra of bounded operators on some finite dimensional Hilbert space B(C d ), 
also known as d-by-d complex matrices. This is also a vector space, which we can equip with 
the Hilbert-Schmidt inner product 

(A,B) := Ty(A^B). 

We have used this scalar product before, so as to combine a density matrix and an effect 
to obtain a probability. Indeed, the density matrix p defines a state which associates the 
probability 

p = Tr (pE) = ( p , E) 

to the truth value of the proposition defined by the effect E. 


Now consider two Hilbert space that of Alice: Ba and that of Bob: Bb- A channel from Alice 
to Bob: £ : B(Ba) —> B(Bb) is in particular a linear map and hence has an adjoint £f with 
respect to the Hilbert-Schmidt inner products on both algebras: 

{£(p),E) = (p, £*(£)} 


or, explicitely, 


Ti(£(p)E)=1r(f>£'(E)). 


The quantity on the left hand side of those equations has a physical intepretation: it is the 
probability that the state £(p) of Bob associates to the effect E of Bob. But we see here that 
is also has a different interpretation from the point of view of Alice: it is the probability that 
Alice’s state p assigns to her effect £\E). 


It is important to note that the adjoint map £' exchanges image and preimage, i.e. £f : 

B(B b ) —► B(B A ). 


In this document, we refer to £ as the channel in the Schrodinger picture, which acts on states, 
and to £t as channel in the Heisenberg picture, as it acts on effects. 

A familiar example is given by unitary channels. If £(p) = UpU\ then £' (E) = U' EU . Indeed, 
Tr (E£(p)) = Tr (EUprf) = Tr(U^EUp) = Tr (£*(E)p). 


We have already established the properties that the map £ needs to be a channel, but what 
do those properties correspond to in terms of £'! One could simply derive them directly from 
those of £. But, equivalently, one may also follow the same principle used to derived the 
properties of £, and demand that £1 preserves convex combinations of effects, which have the 
same ignorance interpretation as for convex combinations of states, and that it maps effects 




3.2 Heisenberg Picture 


33 


to effects, even when acting on part of a larger system. The first condition implies that be 
linear. Since effects are defined in part by the fact that they are positive, the second condition 
implie that be completely positive. 

Hence, both £ and £^ are linear completely positive maps. The only difference is that whereas 
£ must be trace-preserving, £ { must be unital, i.e., £\t) = 1. 

Indeed, for all operators X, we have 

Tr(£ f (l)X) = Tr(l£(X)) = Tr(£(X)) = Tr(X) = Tr(lX). 

Since this holds for all X , this implies 


£\l) = 1 . 


This leads us to a more general definition of channel which holds in more general terms (in¬ 
cluding in the abstract framework of C*-algebras): 


Definition 3.4 (Channel, more general). 

A channel transferring information from the system defined by the algebra A\ to that defined 
by the algebra Ai is represented by a linear, completely positive map 

£^ \ A2 —y A\ ■ 

which is also unital: £t(l) = 1, and meant to be applied to effects. 


In order to avoid confusion, we will typically use the symbol f when using the Heisenberg 
picture, even if there is not corresponding Schrodinger representation. 

Again, it is important to note that in the Heisenberg representation, channels “run backward” 
in time. A consequence of this is that they must be composed in the opposite order from the 
Schrodinger representation, since 

(£ l0 £ 2 y = £\o£l 


The reason this representation is preferred is that for a general C*-Algebras A, effects are 
elements of A, but states are just linear functionals on A. In fact, states can be seen themselves 
as channels from the system defined by the one-dimensional algebra C to that defined by A, 
which, in the Heisenberg picture, is given by a linear completely positive unital map <p : A —> C. 
(This hints at the fact that we are secretly working in the category whose objects are C*- 
algebras, and whose morphisms are unital CP maps). In this language, the effect of & : A 2 —> 
Ai on the state (p : A\ —> C of the system A\ is given by composition: p> is mapped to the new 
state <p' : A 2 —■> C defined by 

tp' — tp o . 
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Figure 3.1: An observable is a transfer of information from the system being observed to a 
pointer recording the result of a measurement. 

3.3 Observables as channels with classical output 

A type of channel which will play an important role in this lecture are those representing a 
transfer of information from any system to a classical system. 

For instance, such channels naturally represent generalised observables of a system. Indeed, 
in the most general terms, an observation is a process by which information about the system 
being observed (represented by the algebra A) is transferred to a measurement apparatus which, 
crucially, is a classical system. Let hi be a set representing the relevant degrees of freedom of 
the classical measurement apparatus (say the possible positions of its “dial”). 

Lienee the following definition: 


Definition 3.5 (Observables). 

An observable of a theory defined by the algebra A is any channel from that system to a classical 
system with pure states Q, represented in the ffeisenberg picture by a unital CP map 

A f : L( Q) ->• A. 

The set O represents the possible values of that observable. 


If this map has a Schrodinger representation 

A :A->L(n), 

then it maps a state p of A to a probability distribution A (p) over the possible set of values O. 


3.3.1 Example: Observable associated with an effect 


Every quantum effect E is naturally associated with an observable: the channel represented in 
the Schrodinger picture by the TPCP map X E : B(T-L) —> L(Q), with O = {false, true} and 

A e(p) = Tr (Ep) |true)(true| + Tr((l — E)p) (false)(false), 
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where we used the orthogonal vectors |true) and | false) to represent the two pure states of the 
classical apparatus. In the Heisenberg picture, it maps classical effects / : O —>■ [0,1] to 

*!(/) = /(true) E + /(false)(1 - E). 


When the system is quantum: A = BfH), we call a classical-valued channel a quantum-to- 
classical channel, or QC-channel. As an example we will determine the general structure of 
such QC-channels in the case where % is finite dimensional and O finite. This allows us to 
represent it as the CPTP map 

A : Bin) -G L(fi), 

From linearity of X we can infer: 


x (p) = {fi(p)}isn = Y^fi(p)\i){i I with f : B{U) ^ C. 

i 


Theorem 2.2 tells us that there exists A* G B{n) such that fi(p) = Tr (pAf). The positivity of 
the map X requires that for all p. Tr(Ajp) > 0 which implies that Aj > 0 for all i. Moreover, 
the trace-preserving property of X says that for all operator A", 


Tr (X(Y)) = ^Tr(AW) = Tr(F) = Tr(lT'). 

i 


Using in particular Y = \i)(j\, this equivalently implies 


E^ = i. 


Hence, the fact that the QC-channel X is linear, trace-preserving and positive implies that it 
is represented by a set of operators A* > 0, i G H, such that A = 1 as 

X{p) = {Tr (Ap)}.^ = E Tr (A/#)(i|- (3.71) 

l 

These operators Aj are in fact effects of the quantum system. Indeed, they automatically satisfy 
0 < A* < 1, and are used to define the probabilities Tr(pAj) which make up the classical state 

X(p). 

Let us show that the complete positivity of X does not impose any extra condition on these 
operators. Any operator W of the extended algebra BfH) ® BfH') can be written as W = 
j A j ® Yj. The action of X ® id on it is 

= E r ^(^i)l*)(*1 ®Yj = E<fc|V / AA J y / A{|A;)|i)(i| ® Yj 

ij ij k 

= T,\i)(k\\fA i X j jA i \k)(i\®Y j 

ijk 

= E (l*><fc| A ® >)(E Xj ® Yj) (f(A t )\k){i\ ® 1) 

ik s " j " --- ' 
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Therefore, we find that X ® id can be represented as (X ® id)(W) = E ik WE} k , which is 

manifestly positive when W is positive. 


In summary, we find that a quantum-to-classical channel between finite systems X : BfH) —y 
L(Q) is defined, according to Equ. 3.71 by a set of quantum effects Ai, i £ 0 such that J2i A = 1- 
A set of effects fulfilling this property is also a special instance of a positive operator-valued 
measure (POVM), a terminology which will become more meaningful when we allow the target 
classical system to have continuously many pure states. 


Consider for example the standard notion of an observable as hcrmitian operator A = Af The 
spectral decomposition theorem tells us that we can write A = ct l P l . Instead of labelling 
the spectral projectors Pi by an integer i, we could also directly use the distinct eigenvalues cq 
themselves as labels, and write 

A = aP a , 

where 12 = {an,..., a at} and P ai = Pi. If the state is p, measuring the observable defined by A 
puts the classical measurement apparatus in state a € 12 with probability 

p a = Tr (pP a ). 

Hence, the observable defines a map from the quantum state p to a probabilistic state of the 
measurement apparatus with pure states 12. This is the QC-channel Xa '■ BfH) —> L(12) with 

X A (p) = {l*(Wy} = E Tr(pfy|a>W, 

aef2 

where the orthogonal vectors |a) are used to label the classical states of the measurement 
apparatus. The set of effects defining this observables are simply the spectral projectors P a , 
and the eigenvalues correspond to the classical pure states 12. 


We see that the old notion of observable defined by a self-adjoint operator corresponds only to 
a special type of QC-channel, or POVM, namely one whose effects are orthogonal projectors. 
We say that it is a sharp , or projective observable. 


Definition 3.6. 

We say that an observable represented by a QC-channel X' is sharp or projective if it maps all 
sharp classical effects E 2 = E to sharp quantum effects: X\E) 2 = X^(E). 


Note that a complete set of projector are automatically mutually orthogonal, which is why 
orthogonality is not part of this definition: 


Lemma 3.1. If Pi are projectors such that Pi < 1, then PiPj = 0 for all i j. 


Proof. For the proof, we just need to show that that it holds for two projectors P and Q such 
that P + Q < 1. Multiplying the last inequality by P from both sides we obtain P + PQP < P 
which implies PQP < 0. But PQP > 0 by construction, so PQP = 0. This implies that for all 
states |-0), (f>\PQQP\f>) = 0. Therefore, HQPIAII = 0 for all |^), which implies QP = 0. □ 






3.3 Observables as channels with classical output 


37 


Despite this fact, the above concept of sharp observable is already more general than that 
associated with self-adjoint operators, in a way which may seem cosmetic at first, but which 
is important to realize. According to our definition, the set of value D can be any measurable 
set, such as, for instance D = M n , something which would require n commuting self-adjoint 
operators to describe. 


3.3.2 Coarse-graining of observables 


There is a natural way in which an observable can be “weakened”, or made less sharp. Consider 
an observable represented by the QC channel X : B(l-L) —> L(£l) (Schrodiner picture). One 
may perform a measurement of this observable on a state p, obtaining the classical probability 
distribution X(p), and then discard information by “scrambling” the output with a classical 
channel n : L(fl) —> L(fT). Effectively, this is equivalent to having measured the new observable 

n o X : B{H) ->L(D'), 

where o denotes the composition of maps: (tt o X)(p) = n(X(p)). We may say that this new 
observable is a coarse-graining of X . It extracts less information about the system than X, for 
any reasonable concept of information. 

For instance, suppose A is a sharp observable: X(p) = Tr(pP Q ,)|a)(a|, where P a are projec¬ 

tors, and we represent 7r by the stochastic matrix 7r a p as n ({p^joen) — {So ^ pala} g&ci'■ Then 
the coarse-grained observable n o X is not in general a sharp observable anymore: 

(tt o X)(p) = {Tr {pEp)} peQI with Ep = 


3.3.3 Accessible observables 

In order to gain some practice with these notions, imagine the following situation: we have 
two quantum systems A and B represented by the Hilbert spaces 1-La and He- There is some 
physical process carrying information from A to B, represented by a CPTP map £ : B(1-La) — >• 
B(1-Lb)- We would like to be able to gain information about system A, however we are just able 
to make measurements on system B. Suppose that we make a measurement represented by the 
QC map X : B(1~Lb ) —<> L(Q). By making this measurement on system B, we are effectively 
measuring the observable 

Xo£ : 

on system A. 

Hence, measuring X on B is perfectly equivalent, in terms of the classical information collected, 
to measuring X o £ on system A directly. This is just the Heisenberg picture: the channel 
represented in the Schrodinger picture by the CPTP map £ sends the observable X to X o £. 
Therefore, we ought to be able to write this action directly in terms of the unital CP map £ { . 
Indeed, this is how the effects defining X transform: if X(p) = {Tr(X Q p)}Q, e Q, we have 

(Xo£) M = {1KJW0,))} = {lK£t(X„))} 
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Hence, £ { maps the POVM elements of X to those of X o £. 

In the sense just explained, we say that the channel £ limits our ability to observe system A, 
and hence the information that we have about A. Specifically, it acts as a constraint of the 
observables of A, since it only allows us to measure effects of the form £ Jf (X), where X is any 
effect of system B. 

As an example consider the channel £ : B(T-La) — > B(Tds ) with 

£(p) = ^TrfpP,)l*>(i|, ( 3 . 72 ) 

l 

where PiPj = SijPi- The adjoint map can be found by considering 

Tr (X£(p)) = Y, Tr(pP f )Tr(X|i)(i|) = £ Tr(p(i|X|i)P) = Tr(p£*(X)), 

i i 

which yields 

£ , (x) = E('lx|i)n. 

i 

Therefore, if we only have access to system B , the only effects of system A that we can indirectly 
observe are of the form J2i x iPi, where Xi G [0,1], which we recognise as the effects of the 
commutative algebra alg(P 1; ..., P n ), which is isomorphic to the classical algebra L({1,..., N}). 
Therefore, from the point of view of an observer of system B, system A looks purely classical! 
This specific channel £ is a simple prototype of the phenomenon of decoherence by which a 
quantum system appears classical. 


3.4 Stinespring dilation 


In the previous sections we have defined the concept of channel abstractly, as the most general 
type of map compatible with the state/effect structure. But we have not explained how such 
maps can occur in practical situations. In this section, we show that channels naturally occur 
in the description of the dynamics of open quantum systems. More importantly, we show that 
any channel has an appropriately unique such representation, which allows for a very powerful 
characterisation of quantum channels which has many consequences. The most general form 
of this result, the Stinespring dilation theorem, is arguably the most important mathematical 
result in quantum information theory. 


We start by explaining how a simple decoherence channel such as the one mentioned in Sec¬ 
tion 3.3.3 can occur in a physically realistic situation. We consider two systems: the “system 
of interest”, S, a two-dimensional quantum system with Hilbert space Ids = C 2 , and an “en¬ 
vironment” E which is also quantum and two-dimensional, with Hilbert space Be = C 2 . We 
imagine that they interact for some fixed amount of time, i.e., that they evolve under a common 
Hamiltonian. The result of this evolution is represented by a unitary map 

U : Ids ® Be —> Ids ® Be- 

By choosing the interaction Hamiltonian, we can obtain any arbitrary unitary map U. For this 
example we consider 

f/ = |o)(o| s ®i E + |i)(i| s ®(;j) E . 
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This unitary is also called the controlled NOT gate in the context of quantum computation, 
where it serves as a basic building block to construct more complex unitaries (See Ref. [T]). 

Here the system S can be thought of as a control bit, since if it starts in state 10), then nothing 
happen to E, whereas, if it starts in state 11), the unitary (? o) is applied to E, which swaps 
states |0) and |1) (a NOT operation in classical logic). Of course, this interpretation of the 
map makes sense only in this particular basis. For instance, in the basis |+) oc |0) + |1) and 
|—) oc |0) — 11), the role of control and target systems are exchanged. 


But let’s keep with the original basis, and ask what happens if it is not the control qubit S 
which starts in state 10), but the target E instead. Fixing the input state of one of the system 
amounts to considering the linear map 


V -.Us ->• E s ® U E 

W U{\^) 0 |0» 


Explicitely, we find 


V = |0)<0| s 0 |0 ) E + |l)(l|s 8 \1)e = |00)(0| + |11)(1|. 

Such a map is a called an isometry , in that, like a unitary map it satisfies V^V = 1, but it 
is not invertible, and VV^ is not the identity (but it is always a projector). It simply embeds 
a small Hilbert space into a bigger one while preserving the orthogonality of vectors. In this 
example; it represents Es as the subspace of Es 0 Ee spanned by 100) and 111). 


Moreover, if the initial state of S is either |0) or jl), the map V simply makes a copy of the 
classical bit encoding this information. Again, this interpretation does not work in a different 
basis, but it indicates that some information is being transmitted from S to E. 


Finally, suppose that we only care about system S after this interaction, i.e., we promise to 
only make further measurements on system S only. In other words we discard E. In the density 
matrix formalism, this amounts to looking at the reduced state of system S in the state after 
the interaction. If the original state of S was represented by p, then after the interaction, the 
state of the joint system is VpV^. Discarding system E yields the final state Tr E VpV^ of the 
system S. 


Altogether, this amounts to the map 

p i y Tt e VpV\ 

which is a channel from S to S. Explicitely, 

Tr E VpV f = ^((*|<g>l)Ep A R t (N)®l) 

i 

= 53((*| ® l)(|i)01 ® \j))p(\k)(k\ ® (k\)(\i) 0 1) 

ijk 

= Y,( i \j)( k \ i )\j)(j\p\ k )( k \ = = E Tr (pKX*l)l*)<* 

ijk i i 


which is of the same form as the channel defined in Equ.(3.72). 
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We see that the unitary evolution transforms a “natural” type of limitation (access only to 
a subsystem) to that of having access only to a commutative subalgebra of observables, as 
explained in Section 3.3.3 for this type of channel. The fact that only classical information is 
left in the system after this interaction is crucially linked to the fact that information about a 
basis was copied between S and E by V, as will be shown in a later part of this lecture. 


In fact, a central result in quantum information theory is that any channel can occur in this 
way. This is stated by the following theorem, which have have here specialised to the finite¬ 
dimensional case in order to make the proof more accessible: 


Theorem 3.1 (Stinespring). A linear map £ : B(Ha) —» B(H E ) is completely positive if and 
only if there exists a Hilbert space He (environment) and a linear map V : Ha —>■ Hb ®He 
such that, for all state p, 

£(p) = Tr E (VpV t), 


or, equivalently, 


£\X ) = V\X®t E )V. 


The map £ is also trace-preserving, and hence a channel, if and only ifVW = 1. 


Proof. This theorem in fact holds when the Hilbert spaces Ha and H E are infinite-dimensional, 
but for simplicity we only prove it here in the finite-dimensional case, where we present a proof 
adapted from Ref. [3]. We immediately see that any map of the form £{p) = Tr E {VpV^) is 
completely positive because it is the composition of two completely positive operations, namely 
p H- V pW and Tr^. 


Conversely, suppose £ is completely positive. Let us define the matrix 

X £ (£ 0 id)(|0)(0|) with |Q) = - |i) 0 |i), 

i 

where |z) forms a basis of Ha- Clearly, X? > 0 since £ is CP and |f2)(f2| > 0. This object X £ is 
commonly referred to as the Choi matrix of £. The good thing about the Choi matrix is that 
we can reconstruct the CP map £ from it as follows: 

£{p) = (li 0 <n| 23 )(£i 0 id 2 0 id 3 )(|n)(n|i 2 ® p 3 )(li ® 1^)23) , . 

= (li ® (n| 23 )(x£: (8 p 3 )(li ® 1^)23) 


where we are now using three different copies of Ha which we labelled by the subscripts 1, 2, 3 
for clarity. This equation can be verified by direct expansion of the definition of the maximally 
entangled state Q). This invertible map between CP maps and positive matrices is known as 
the Choi- Jamiolkowski isomorphism. The particular representation of the CP map that we are 
looking for simply amounts to diagonalising the Choi matrix. Indeed, since X £ > 0, we can 
write X £ in diagonal form as X £ = Yffi= 1 where cq > 0 and | 'tfj ) = S iy Substituting 

this diagonal form in Equ. 3.73 yields the expression 


£(p) = (1 0 (fi|)C\:£ 0 p)(l 0 |fi)) 

= ® ( k \)\A)(k\p\j){A\(^ 0 li)) = X) EipE J, 

ijk i 
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where 

Ei = J2 ® ( k \)\A)( k \ 

k 

are called Kraus operators for the CP map. They are directly related to the operator V that 
we are looking for. Indeed, if we now take to be a hilbert space of dimension equal to the 
rank d of Xg, and {|i)}f =1 any orthonormal basis of He- Then if we define 


V :=Y^Ei® | i) E . 

i 


We obtain 


Tv E VpV' = '£ E -pE! = £(p), 

i 

which is the form that we were looking for. Moreover, if £ is also trace-preserving, then £l is 
unital, which means that £t(l) = V^V — 1. 0 


The isometry V from Theorem |3.1[ together with the Hilbert space He, is called a Stinespring 
dilation of £. An important aspect of this dilation is that it is unique up to a partial isometry 
on He- A partial isometry W preserves the scalar product on its range, but can send some 
vectors to zero. It is characterised by the fact that W'W is a projector, which also implies that 
WW f is a projector. 


Theorem 3.2 (Stinespring, uniqueness). Let He, and He 2 and V± : Ha —» Hb 0 H El , V 2 : 
Ha —t H E ®H E2 be such that 

Vl{X 0 1 £l )W = Vj(X 0 1 E2 )V 2 

for all X G B{H E )- Then there is a partial isometry W : H El —t He 2 , suck that 

V 2 = 


Proof. Consider the span /C of all vectors of the form (A" 0 1 e 1 )I / i|' 0), for arbitrary state \if) 
and operator A". We consider a basis (X i ®\ El )Vi\iff) for i — 1 ,.. . ,n and define W as mapping 
(A f 0 l E f)V\ \f>i) to (X l 0 t E2 )V 2 \'iff). If \<f>) belongs to the orthogonal complement of /C, we 
define W\<f>) = 0. It follows from this definition that W{X 0 t)V\ — W{X® t)V 2 for all X, 
and in particular, 

WV i = V 2 . 

The operator IT is a partial isometry because, due to the fact that both dilations represent the 
same channel, 

(^|T 2 t (A: t 0 i)(y 0 i)v 2 \rf) = (i/>\ v}(x*y 0 t)v 2 \if) 

= (if\Vi (x^y 0 i)w| Tjj) = {ii>\vf(x* 0 i)(y 01)Vi|-0) 

for all X, Y and all \if). Moreover, for all X , Y, 

w(x 0 i)(y 0 i)v x = {x 0 i)(y 0 i)v 2 = (x 0 i)w(y 0 i)vi, 
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which shows that W(X<g> t)\4>) = (X 0 t)W\ft) for all \ft) G 1C. Moreover, if \ft) belongs to the 
orthogonal complement /(V, we have W\ft) = 0 by definition of W, but also, since W' maps 
onto /C, then for all \ft) and all X, (ip\ W(X 0 1)| ft) = 0, which implies W(X 0 1 = 0. 

Therefore, W(X 0 1)|0') = (X 0 1)W\4>') also for all ft G X ± . We conclude that 

W(X 0 1) = (X 0 1)W 


for all X , which implies 


W = 1 0 W 


for some operator W which inherits its isometric property from W. 


□ 
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